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Abstract 

Dynamical breaking of chiral symmetry in effective models of QCD is studied. In- 
troducing a cut-off function or a non-local interaction, the Noether current must be 
modified and thus the Ward-Takahashi identity and the PCAC relation are modified 
accordingly. We point out that the pion decay constant must be defined consistently with 
the Noether current so that the low-energy relations are satisfied. We define the proxy of 
the Noether current for general effective models, which is consistent with loop expansion 
of the Cornwall-Jackiw-Tomboulis effective action. A general formula for the pion decay 
constant in terms of the Bethe-Salpeter amplitude is derived. The effective Pagels-Stokar 
formula is proposed which is useful to estimate the decay constant without solving the 
Bethe-Salpeter equation. 
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1 Introduction 



The program to derive the observed properties of hadrons non-perturbatively in quantum chro- 
modynamics(QCD) has been pursued with great intensity but not accomphshed yet. The 
concept of chiral symmetry and its spontaneous breakdown are among the most important as- 
pects of low-energy hadron physics. The spontaneous breakdown of chiral symmetry is believed 
to be responsible for a large part of the low-lying hadron masses as well as for the emergence of 
octet pseudo-scalar mesons as Nambu-Goldstone(NG) bosons. In order to explain the observed 
hadron spectrum, one also needs small, explicitly chiral symmetry breaking terms, namely, the 
flavor dependent current quark mass terms. 

Dynamical chiral symmetry breaking(DCSB) of QCD has been extensively studied in effec- 
tive models in terms of light quarks. The Nambu-Jona-Lasinio(NJL) modelf^, as an example, 
is a superb model, that presents the main idea of DCSB most concisely. It describes the 
spectrum and properties of NG mesons fairly well.0] The NJL model, however, has several 
shortcomings. It may be too simple to represent some of the key properties of QCD. First, it 
does not describe color confiniement and therefore the meson spectrum above 2Mq threshold 
(Mq is the constituent quark mass) is not reproduced. Secondly, the quark interactions are not 
asymptotically free in NJL. This property will be crucial to describe short distance structure 
of mesons. 

There have been various attempts to derive effective quark theories which are consistent 
with the asymptotic behavior and low-energy chiral symmetry breaking. One of such models 
is the improved ladder approximation proposed by Higashijima and Miransky. 0, ^ This model 
is consistent with the asymptotic behavior of the leading order renormalization group analysis 
of QCD.0 It also realizes DCSB, as is seen in the non-zero quark mass function given as 
the solution of the Schwinger-Dyson(SD) equation. Then according to the NG theorem, the 
Bethe-Salpeter(BS) equation gives a massless NG boson with J^'~"' = O^^.0 The numerical 
results for the pion decay constant and the quark condensate {-ip^) are consistent with other 
analyses. 0] On the other hand, it is known that this approximation violates the axial- vector 
Ward-Takahashi identity. ^ Therefore it looks unsuitable to study the low-energy relations in 
QCD. 
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However, there exists a global chiral symmetry in the improved ladder approximation. It 

produces a Noether current associated with the axial transformation, namely axial-vector cur- 

— A" 

rent J^^, which is generally modified from the QCD form ^/'7^75 — ip. Accordingly the pion 
decay constant defined by 

tp,fJa^:={0\J^,m^^ip)) (1) 

must be modified. The Ward-Takahashi identity for the modified axial-vector current is satis- 
fied. 

The improved ladder approximation is considered as a typical example of the effective model 



which contains a non-local interaction term, while various similar models are proposed. [0, 10 



|TI], |T^, [13], |T^, jTB], |T^, [IB], |TP| Among these effective models some require no modification 



of the axial- vector current. However, the corresponding Bethe-Salpeter equations become so 
complicated that practical calculation is rather difficult in such models. In Ref.0, the trun- 
cation in the Gegenbauer polynomial expansion is employed. This procedure may violate the 
low-energy relations. In Ref.^], to preserve the axial- vector Ward-Takahashi identity, the 
running gauge parameter is employed, which makes practical calculation hard for finite quark 
mass. In Ref. [^], the SD and BS equations are solved consistently. But as the wave function 
renormalization constant Z2 deviates from 1, the vector current or the charge operator should 
be modified. 

On the other hand, the improved ladder approximation is simple and useful to treat numer- 
ically. The aim of this paper is to study a systematic approach which repairs the low-energy 
relations in the effective models and makes it possible to study the low-energy relations in 
QCD. We propose to define the axial- vector current consistenly with the approximation in the 
effective model, and describe the low-energy relation in terms of the solutions of the SD and BS 
equations. We derive a formula for the consistent pion decay constant. There have been similar 
studies of this problem in the literature |T^, where the NJL models with smooth cut-off 
regularization are studied. However our approach is more general and systematic. We study 
the SD and BS equations in the loop expansion using the Cornwall-Jackiw-Tomboulis(CJT) 
effective action formulation.]^] We show that the BS equation has a massless bound state 



solution corresponding to the NG boson. We also give a formula for the decay constant which 
is consistent with the rainbow-ladder approximation. 



In Sec.^ we present how the Noether current is modified due to the loop momentum cut-off 
and/or the non-local interaction. The modification of the axial Ward-Takahashi identity is also 
discussed. 

In Sec.^ we employ the CJT effective action to study the consistency of the SD and the BS 
equations. The NG theorem and the pion decay constant are studied in terms of the effective 
action. We further consider the explicit chiral symmetry breaking due to the finite quark mass. 
Three cases in the patterns of the local chiral symmetry breaking of the effective models are 
studied in detail. 

In Sec.^, we present a general formula of the pion decay constant in the loop expansion of 
the CJT acion. The formula is given in terms of the quark full propagator as the solution of 
the SD equation and the pion BS amplitude. We also propose a Pagels-Stokar type formula, 
which gives the pion decay constant in terms only of the quark full propagator, namely, the 
mass function B{q^). 

In Sec.|^, we employ a simple numerical model to check our analytical results. We confirm 
that the modifications of axial-vector current and the pion decay constant are significantly 
large, while our new Pagels-Stokar type formula gives a good approximation. 

A conclusion is given in Sec.^. 



2 Noether current and Low-energy relations 

Before the disscussion about the chiral symmetry and low-energy relations in the approximated 
SD and BS equations, we consider general aspects of chiral symmetry and low-energy relations 
in effective models. In this section, we show how to modify the Noether current and the Ward- 
Takahashi identity from the original form in QCD. This gives a perspective in the following 
discussion for the case of the rainbow-ladder approximated SD and BS equations. 
We consider a general effective model of QCD whose lagrangian density is given by 

:= £free[V^,V^] +/:int[^,?] (2) 
Aree[^,V^ := ^ f {d^){l^ - TTlo)^ (3) 

where the quark field ip is a, column vector in the color, the fiavor and the Dirac space. Here the 
function /((^) oi C, = is introduced as a cut-off regularization function in order to regularize 
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the ultraviolet divergences coming from quark loops. The reason we introduce the cut-off 
function at the lagrangian level is to preserve the consistency of the SD and BS equations. If 
one uses the regularization that is inconsistent between the SD and BS equations, the low- 
energy relations based on the chiral symmetry should be violated by the regularization. The 
function f{() should satisfy f{( = 0) = 1. For ( ^ Mjv^ /(C) should diverge sufficiently fast 
so as to regularize loop integrals. For example, the sharp cut-off is given by 

/(C) = l + M^(A^v-C), M^oo (4) 

and the Gaussian smooth cut-off is given by 

/(C) = l + Mexp(^^^. (5) 

A caveat of this procedure is that a general choice of f{() may cause a difficulty in canonical 
quantization. Since our purpose is not to study the non-local field theory, we simply employ the 
path integral formulation with the action integral regularized by /(C)- Although this treatment 
is not rigorous, it is sufficent in the present disscussion. We thus maintain the consistency 
between the chiral symmetry and the regularization. 



/(C) also determines a scale Auv at which the bare quark mass mo is evaluated. |g, |TT 
To compare the bare mass rriQ the renormalized quark mass in QCD, one must impose 
a suitable renormalization condition. In general mo is a diagonal flavor matrix, i.e., tuq = 
diag(m„, m^i, m^) for Nf = 3. But in this paper we deal only with a flavor independent mass 
or the SU{3) limit. The generalization to the flavor dependent masses can be also done. 
Although there is a difficulty in choosing the center of mass coordinate for the quark-antiquark 
bound states physically, the low-energy relations hold for any choice of the center of mass 
coordinate. 

A general non-local 4-quark interaction is written by 
C,ntbPM^) ■= -If /C'"'"'""'(p,p';g,g') 

A Jpp'qq' 

X ?^(p)^n^'(p')V^n(g)^n'(?')e-^(^+^'+^+'^')^ (6) 

where denotes / and the Fourier transformation of the quark field is defined by 



tfjip) := J dWP"V(a;), (7) 
^{p) := / d^xe'P'^ijix). (8) 



The indices m,n, - ■ ■ are combined indices m := {a,i, /), n := {b,j,g), ■ ■ ■ with Dirac indices 
a,b, - ■ ■ and color indices ■ ■ ■ and flavor indices f,g,-- ■■ We do not consider the quantum 
correction and renormalization for the interaction itself, because we assume that effective models 
have already contained such effects. /C"''"''""'(p,p'; g, g') is an interaction kernel which we can 
assume without loss of generality 

We further assume that the interaction is invariant under the global SU{3)l x SU{3)ii trans- 
formation 

/C'"™'"'(p,p';g,g')(^Y)^n' = (2y)n,/C™"^' g, g') (10) 

/C'"™'"'(p,p';g,g')(^75Y).n' = -(^75 Y)n./C"^'"' '"'(p,p'; g, g') (11) 

where A° is the Gell-Mann matrix [a = 1, ■ ■ ■ , 8) for Np = 3 with normalization condition 
tr(^) [A"A'^] = 25°'^. Thus the violation of the global chiral symmetry comes from the quark mass 
term only. We here concentrate only on the chiral symmetric 4-quark interaction for simplicity. 
The generalization to the multi-quark interaction is straightforward. Some effective models 
in Ref. p4|, |16|, contain the interaction which violates the global U{1)a symmetry, and our 



approach can be generalized to such models as well. 
Under the infinitesimal axial transformation 

^{x) ^ ^'{x) := (1 + Z75y r (12) 

the action 

S[ij,i^] := J d'^xC[ilj,i^]{x) (13) 
changes by 

rf^xr(x)|v^(x)(> fCd'HY + /(^')^75y)V^(a;) 

- 2moz?(x)^^^5^^±^75Y^(x) (14) 
f/C""^' ""'(P,p'; -P-P'- q', q') - /C"^'">"'(p,p'; g, -p - p' ~ q)) 

pp'qq' \ / 
X?^^(p)^r„'(p')?^n(g)(^75Y^)n'(g')e-''(^+^'+''+^'^4. 



One finds that the third term in the right hand side(RHS) of Eq.(14) vanishes if the momentum 
dependence of the interaction kernel is such that 

This relation is satisfied when the kernel is generated by one gluon exchange whose coupling 
depends only on the transfer momentum.]^, [10|, [12], |13[ Furthermore if the cut-off function f{Q 
is identically 1, we obtain the operator identity 

d^J^^ix) = 2moJ^{x) (16) 

for 

Jlix) := V^7^75^^(x), (17) 

J^{x) := V^^75^^(a;). (18) 
For general /C™'"''"'"'(p, p'; g, q') and /(C), Eq.([T^) is written as 

A5^[^,^] = j (fxe''{x)\d''J^^ - 2mo/5"}(x) (19) 



with 



^5,(^) := /5"m(^)-^.°(^)- (20) 

Jp,q 

+ ip- q),{i> - i)f2{-p\ -g'))75^^(g)e-^(^+'')^ (21) 



i{p + p' + q + q')f, 
pp-qq' {p + p' + q + q'y 

X f/C'"'"''""'(p,p'; -p-p' - q', q') - /C"*'"''""'(p,p'; g, -p - p' - q) 



x^^(p)^™'(p')V'„(g)(^75Y^)n'(g')e-^^^+^'+^+'^'^^ (22) 

It{x) := ?(x)^^^^^^±^z75y^(x), (23) 

2 2^ /( V) + fi-l') r, 2 2^ /( V) - /( V) .... 

h{-p ,-q)-= ^ , f2{-p , -g ) := — ^(^2 _ — • (24) 

J^^(x) is the Noether current associated with the axial transformation Eg . ([l^) and is conserved 
in the chiral limit mo 0. For non-zero mg, it satisfies the operator identity 

d^^J^^ix) = 2moI^{x). (25) 



Because the effective lagrangian contains non-standard momentum dependencies, the Noether 
axial-vector current is non-local and cannot be calculated in the standard procedure. To avoid 
confusion, we call Eq.(^^ the naive axial- vector current and call Eq.(^) the true axial- vector 
current. 

When chiral symmetry is broken dynamically, the Schwinger-Dyson equation has a solution 
with non-zero quark mass function and the Bethe-Salpeter equation gives a massless pion state 
which must appear as a pole of the Noether current according to the NG theorem. We define 
the effective decay constant /^r associated with the Noether current by 

iPb.L := {^\JtMP) (26) 
which is compared to the naive decay constant /^r defined by 

iPB.f. := {Q\JtMP) (27) 
where |P) denotes a pion state with normalization condition {P\P') = {27c)^2PboS^{P — P') 



and Pb^i '■= {y + P^, P) denotes the on-shell momentum. The matrix element of Eq.(^ 
between a pion state {P \ and the vacuum |0) gives 



where 

P P 

q_B := q - q+B ■■= q+ ^ (29) 

and the BS amplitude x(9; Pb) and its conjugate are defined by 

xiq;PB) ■■= e'''^'' J d\x-y)e'^^^-y\0\T^ix)i^{y)\P), (30) 
x{q;PB) := e-'^-^|rf^(y-a;)e^«(^--)(P|T^(2/)^(x)|0), (31) 



with X := It should be noted that Eq.(^) is an exact relation representing Eq . (^51) , where 
exact means it holds for any itlq. 

The Ward-Takahashi identity can be derived by the standard path integral procedure 

-td''{0\T^{xmy)J^^iOm = -2zmo(0|rV^(x)^(y)/nO)|0) (32) 
+ i^,^{0\T^P{xmy)\0)S\x) + {0\Tij{x)i^{y)\0)t^,^5\y). 



Note that in Eg. (|32|) we replace the T* product to the T product. (This proceduce is not 
rigorous but is justified when J^^(O) and /^(O) are local currents.) From this identity, one 
obtains the NG solution of the BS amplitude in the limit of mo — > 0, 

x(g;0) = i{z^75,^F(g)} (33) 

where Spiq) is the full propagator of the quark. For a small quark mass mo, one expects that 
Eq . (^31) gives an approximate solution of the BS equation. An approximation of the exact 
relation Eq.(p8D is obtained by substituting Eq. (|33|) into Eq. (|28D 

^ -2mo(V^^)o, (V^^)o := - / f{-q'MSF{qU=o]. (34) 

Jq 

The chiral condensate {ipip)o is considered to be evaluated at the scale Auv- To compare it the 
renormalized one in QCD, one must impose a suitable renormalization condition. This is the 
Gell-Mann-Oakes-Renner mass formula. 

In conclusion, we have formally proved that the low-energy relations are satisfied if we rede- 
fine the axial-vector current and the decay constant f-^ according to the global chiral symmetry 
of the effective lagrangian. In the next section, we show how the above relations are described 
by the quark propagator and the BS amplitude in the framework of the rainbow-ladder ap- 
proximation. Later, we will see that the difference between f^^ and /^r is significantly large in 
the improved ladder model. 



3 Analysis with CJT Action 

In this section, we consider the chiral symmetry property of the SD and BS equations when 
they are truncated. The similar discussion has been already done in the case of QED or for a 
specific interaction kernel which satisfies Eq. (|l^) . 12^, |2l|] However, it is non-trivial in the case 
of the general non-local interactions. Especially it is difficult to find the formula of the decay 
constant f^^ because the Noether current is complicated. For the systematic discussion, we use 
the Cornwall, Jackiw and Tomboulis (CJT) effective action formulation. 

The Cornwall, Jackiw and Tomboulis effective action formulation is one of the most powerful 
and useful methods to derive the SD and BS equations consistently with the chiral symmetry. [^, 



9 



25, 2^ We study the chiral property with this formulation following Munczek. |21] The CJT 
action is given by 

T[Sf] := tTiLniSp] - tTT[So'SF] + Ti.oASf]- (35) 

The last term of Eq. (|35|) , iFioop [Sp] is given by the sum of all Feynman amplitudes of 2-particle 
irreducible vacuum diagrams with two or more loops in which every bare quark propagator 5*0 
is replaced by the full one Sp. We will show that the truncation in the loop expansion preserves 
properties of the chiral symmetry. 

Using the CJT action, the SD equation and the inhomogeneous BS equation are derived by 
ST[Sf] 



6SF{x,y) 



0, (36) 



1 5'^r[>s' ] 

-T^ 7 ri~i\^C-n'm'm"n"iy'^''^ ^' V") = Sm"mSnn"S{x" - x)5{y - y") (37) 

Z ODFmn{X,y)ODFn'm'{y ,X ) 

respectively, where the repeated indices are summed and the repeated arguments are inter- 

(2) 

grated. GQ'j^,^,^„n„{y'x' ; x"y") is the two-body connected Green function defined by 

Gclmm'n'iy^'^^'y') ■= (0 |T^n (Z/) (x) (x')^„/ (?/') 1 0) 

- (O|T^„(y)V^„(x)|O)(O|TV>^.(a;')V^n'(2/0|O). (38) 

(2) 

By representing G(^.„,^,^„„„(?/'x'; x"y") in the spectral form and taking the pion pole term, we 

(2) 

express Gq in terms of the BS amplitude of the pion as 
n(^) (..^. ^^/^ _ f p-i{P-PB)x+i(p-PB)x' 

'^C;nmrn'n'\yX, ^ y ) — I t, 

^^Xurnhj^X^^ ^ i?_,„,(y^; ^V) (39) 



with 



X:=^, r:=^, (40) 



P=(Po,P), Pb = WM^ + P',P). (41) 



The regular term Rnmm'n'iyx; x'y') denotes the contributions from excited states. The BS 
amplitude Xnm{y, x; Pb) in Eq. ([39|) is a solution of the homogeneous BS equation, given by 

'^^ ^Xn'm'(y',x';PB) = 0. (42) 



5SFmnix,y)5SFn'm'iy', x') 
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In order to obtain the normalized BS amplitude, we go back to Eq.(|3Y|) and determine the 
normalization. 

The question is whether the chiral properties of the effective model are preserved when var- 
ious approximations are taken into account. To answer this let us consider a local infinitesimal 
axial transformation of SF{x,y) 

Spix,y) -> S'Ax,y) := (1 + ^75 Yr(x))5^(a;, + ^75 y^'ly)) (43) 

Under this local transformation, the change of the CJT action is 

A5r[^^] = ,/^^^^| ., {^75^^^ S^U^iy', x') (44) 

with 

{^75 SF}{y,x) := ^^,^e^{y)SF{y,x) + SF{y,x)t^,^e^{x). (45) 
Then the following equation holds: 
.(2) ....... 



^(2) / // // N 5^r[5'i?] na o ^ II / 

= Gb'm"n"nm{x V '^V^Jj^ / , xcy TZTZa^^^^^^ ,SF}m'n'{x,y 

Oopnm [y,X)0 OFm'n' {X ,y ) ^ 



+ ^/^^t! ,M S-,rn''n>'Ux^'v'-. ) (^5 (y') ).n' 

OOFra'n'KX ,y ) I ^ 

+ (^75yr(x'))™.4;L"n'.n'(^"/; ^'v')] 

= ^{^75Y^",5^}™"n"(a;",2/")- (46) 

Here for the last equality we use the SD equation (^) and the inhomogeneous BS equation 
iW^). Eq . (^61) is a key equation in our study of the system of the SD and BS equations. 

In the following, we consider three cases in which the local chiral invariance is broken in 
different ways. First, we discuss the naive case where neither cut-off nor momentum dependent 
interactions are included. Second, we take into account the loop cut-off function. Finally we 
consider the case where the interaction kernel is not locally chiral invariant. In each case, we 
will show that the "pion" becomes massless in the chiral limit, and present the formula for the 
pion decay constant. We further show the PCAC relation in each case. 
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3.1 Naive Case 

First we consider the case that the cut-off function is not introduced and the interaction does 
not modify the axial- vector current. In this case, truncation to an arbitrary subset of diagrams 
in Fioopi^F] does not violate the local chiral invariance.||27|| Therefore, the following results are 



valid in each order of the loop expansion of rioop[5'i?]. The first term zTrLn[S'ir] in Eg . (1351) is 
also invariant because 

A,iTiLn[SF] = zTr[5^H^75 Sp}] = 2«tr[«75y] J t/^x^ (x) = 0. (47) 

Then only the second term — ■iTr[5'o"^5'i7'] in Eq.(p5D contributes in the left hand side(LHS) of 
Eg. (^61) , and gives in momentum space 

= <n45HP-^)).„„„„ + <s,(p + |)n5y)„„„„. (48) 

Here the Fourier transformation of G(j.^i^,^^{x'y'] yx) is defined by 



A" 



JpqP 

Using Eg.(|9D, Eg.(|8|) becomes 

:Xm"n"{p; i^B)Xnm(9; Pb) p . p. 

+ ^m"n"nm[P, Q, -T J 



with 

Xmn{x,y;PB) = e-^^-^ /e-^''(^-^)xmn(g;PB), (51) 

Jq 

Xnm{y,x;PB) = e^''^^ /e-*^(^-^)x„^(g;PB). (52) 

Jq 

Chiral Limit 

In the chiral limit rriQ —>■ 0, the CJT action is chiral invariant. We confirm that the "pion" 
becomes massless in this limit because if M,r 7^ LHS of Eg.(^) diverges in the limit ^ M^, 
while RHS converges. In this case LHS of Eg.(|50|) becomes 

LHS of Eg.(50) = J ^tT[x{q; PB)tl5Y^]Xm"n"{p; Pb) 
r A" 

~l~ / Rm"n"nm (p,g;P)(z75— f)„.. (53) 

Jq Z 
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If we take the soft limit P after taking the on-shell hmit P Pb (i.e. Pq ~^ v -P"), the 
second term of Eq. ( |53D vanishes and the first term becomes 

1st term of Eq.(53) if7TXm"n"{p; 0), (54) 

with 

U = \im^ ^ tr[x(g; Pb)^75 ]• (55) 

On the other hand RHS of Eq. (^) becomes 

RHS of Eq.(50) -> ^{^75 y, ^F(p)}m"„" (56) 

in the same hmit. Then we obtain the equation 
1 A° 

Xmniq; 0) = -p{«751^, SF{q)}mn- (57) 

This is the same equation as Eq. (|33| ) which was derived from the axial- vector Ward-Takahashi 
identity. It should be noted that Eq . (^Tf) is valid in the truncated SD and BS equations in 
which Fioopi^i?] is expanded to a finite number of loops. 

Finite Quark Mass 

When mo > 0, the first term in the brackets of Eq.(pOD will diverge in the on-shell limit P Pb- 
Then the equation 

|^tr[x(g;PB)«75Y(2mo + fB)] =0 (58) 

must be satisfied and we obtain 

M^U = -2m^i jii[x{q; Pb)i^^] (59) 

where /,r is given by Eq. (pSj) . This equation is again same as Eq. (pSD . Using the approximation 
Eq. (p7D , we obtain the GMOR relation 



Mlfl ^ -2mo{^^)o, (^V^)o := - / tT[SF{qU=o]. (60) 

Jq 

Note again that any truncation of rioop[5'ir] is guaranteed to satisfy all the above equations. 
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3.2 The Case of the Cut-ofF Regular izat ion 

In order to regularize loop integrals, we introduce the cut-off function in the kinetic term 
of the lagrangian as in Eq.(^). As far as this regularization is taken, the interaction term 
of the effective action is not modified and therefore the chiral invariance of rioop[>S'_F] is not 
violated. When the cut-off function is introduced, the same procedure can be apphed with the 
replacement 

275y (2mo + f ) (61) 

. A"/ /(-d) + /(-g^) 2w 2w \ 

P P 

q^:=q-—, q^:=q+—. (63) 



As a result the decay constant is given by 

A"r/(-g^) + /(-g^) 



L = ^ tr[x(g; PB)^7,Y{ 2 + ^^^'"'^^ " ^^'^-^ 



(64) 



This coincides with the definition (pGl). Therefore the exact relation Eg. ([59|) is modified to 



M'jl = -2m,^ / ^^^^^^^-^^^tr[x(g; Pb)i.^] (65) 



and the definition of quark condensate is modified to 

(^^)o = - / f{-q^MSF{qU,=o]. (66) 

Jq 

Under these modifications, Eqs. (|57|) and ( |60| ) hold. 

3.3 The Case of the Non-local Interaction 

Here we consider the case that the interaction modifies the axial- vector current. To simplify 
the argument, we omit the cut-off function f{Q in this subsection. We employ the two-loop 
approximation of Fioopi^'i;'] such that 

Tioopl^i^] = -^ld^x]C"'^''^^^^^''^td^„td,,;tdy„tdy,) (67) 

X [SFm2mi{x2, Xl)SFn2niiy2,yi) — S Fm2niix2, Vl) SFn2mi{y2, Xl)] \^ 
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where the symbol * means to take xi,X2,yi,y2 — x after all the derivatives are operated. Note 
that this approximation corresponds to the rainbow approximation in the SD equation and the 
ladder approximation in the BS equation. Then Eq.(H6D reduces in momentum space to 



A" 



G'cU"n"nm{P, P) (^5 Y(2mo + f ) + E^{q; P)] 
P P 1 A" 



ln,n'm' 



/c 



I J jQn'm',mn 



I J j^n' n,mm' 



P P 

-k + P,k;-q- -,q- -) 

P P ^ A" 
'i-k,k + P;-q-—,q- — )|(i75 Y^F)m'n'(/i;) 

P P 
-k,q- -■,-q- j,k + P) 



,mm 



P P 

-k + P,q-—;-q-—,k) 



A" 

^75^(27710 



Substituting Eg . (|39D , Eg . (|68|) becomes 

• Xm"n" X,nm 



P^-Ml + 



le 



~l~ Rm"n"nm 

{p, q; P) 



(68) 



(69) 



+ E^{q- P) 



(70) 



2 ^ 2 V ' 2 2 

It is easy to see that E^lq; P) becomes zero in the limit P ^ 0. This is expected because LHS 
of Eq.(|68D or Eq.(^) should be finite in the chiral limit as RHS is. It is also easy to check that 
E"(g; P) vanishes if /C"^™''""'(p,p'; q, q') satisfies Eq.(|l5|). 



Chiral Limit 

In the chiral limit, the global chiral symmetry breaking implies the existence of the NG boson 
again. Then we obtain 

LHS of Eq.(70) = ^ -^tr[x(g; Pb)(^75 yf + E^{q- P))]x„^"n"(p; Pb) 
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+ / Rm"n"nm{p,q]P){il5-^f)mn- (71) 

Taking the soft limit P — ^ after the on-shell hmit P ^ Pb (i.e. Pq ~^ Vp^) the second term 
of Eq. ([TTl) vanishes and the first term becomes 

1st term of Eq.(71) ifT,Xra"n"{p] 0) (72) 

where we use a (new) definition of /^r, 

h := ;nn^ /^tr[x(g; Pi.)(^75 yf + i?°(g; P))]- (73) 

On the other hand RHS of Eq. ( |7D| ) becomes 

RHS of Eq.(70) ^ ^{^75 y, ^F(p)}m"n" (74) 

in the same hmit. Thus we obtain the relation 
1 A" 

Xmn{q\ 0) = ^{^751^, SF{q)}mn- (75) 

This is the same equation as Eq. (|33|) which was derived from the axial- vector Ward-Takahashi 
identity. Therefore we conclude that the definition of f.^ in Eq. (|73D is equivalent to the definition 
Eq.( 



Finite Quark Mass 

When mo > 0, the first term in the brackets of Eq.(^) diverges in the on-shell limit P Pb- 
Then the equation 

1^ tr[x(g; PJ,){^^^^{2m, + f b) + E"{q- Pb))] = (76) 
must hold. From Eq. (|73|) we obtain 

r2 " - A 



M'J^ = -2mo2 y^tr[x(g; Pb)!,,^]- (77) 

This equation again coincides with Eq.(P5D when f[Q = l and implies that our choice of /tt in 
Eq . ([73D is consistent with Eq . (p6D . The GMOR relation is written as 



M^V. ^ -2mo(V^V^)o, (^V^)o := - / tr[Si.(g)„„=o]. (78) 



9 
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Thus we have proved that the truncation of FioopI'S'f] preserves the low energy property of the 
effective model if one uses the appropriate formula of the decay constant given by Eq.(|73|). 

Up to now, we neglect the cut-off function. Introducing both the effects of the cut-off 
function and the non-local interaction, the formula for the decay constant in the approximation 
Eq.(p7D is given by 



(79) 

It should be noted here that we have defined in Eq.([7^ so as to reproduce the low-energy 
relations Eqs.(^^ and ([77D. The definition Eq. ([75D is obtained from Eq.([2B|) directly by taking 
the following approximations 

(0|T^V^|P) ^ (0|T^^|P)ladder, (80) 
(0|Wt^^|P) ^ E(0|^^V^|0)rainbow(0|TV;?^|P)ladder. (81) 

4 General Formula of the Decay Constant 

So far we have treated the lowest order (rainbow-ladder) approximation in the CJT action. 
When we proceed to the higher order terms, the two- loop formula Eq.(^) should be modified 
to an appropriate form. Then we can prove that the truncation of Fioopi^'F] preserves the low 
energy property of the effective model. We note that the consistency is guaranteed in the loop 
expansion of the CJT action formulation. Further approximations inconsistent with the loop 
expansion will violate the low energy properties. An example is to take the leading terms of 
the Chebychev polynomial expansion of the BS amplitude, although the numerical result shows 
that the violation of the PC AC relation is generally small. |T2|, [TB| 



In order to derive a general formula of the decay constant, we consider the transformation 
property of the CJT action. Under the infinitesimal local axial transformation, the change of 
the classical action can be written as 

A,S[^,^] = [ d'xe"{x){d^J^^[ijM^) - M'^[ij,^]{x)). (82) 
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Here M"'[iIj,iIj]{x) comes from the globally variant terms, such as the quark mass term and 
J^^[ip , 'ip]{x) is an effective axial- vector current. 

The change of the CJT action corresponding to the classical action must be written as 

A,T[Sf] = I d''xe''ix){d^J,-^[SF]ix) - M''[Sf]{x)). (83) 

Again A4"'[Sf]{x) comes from the globally variant terms and J'^^[Sf]{x) is a proxy of the 
effective axial-vector current in the CJT action. 

When the effective axial- vector current J^^['ip,ip]{x) has a non-local interaction term, the 
exact proxy J'^fj_[S f]{x) is an infinite sum of the Feynman amplitudes which come from the 
expansion of rioop[>S'F]- If one truncates the expansion of rioop[5'_F], the approximated proxy 
J'51[Sf]{x) is built of a finite sum of the Feynman amplitudes. 

For example, in the lowest loop approximation Eq.(|67D and with the cut-off function f{() 
we obtain 

d^J,l[SF]{z) = J rf^x^e-^^(^-^)/(V){tr[|^275Y5^(^,a;)] -tr[|^5^(a;,^)z75Y]} 

( A" A" 1 

r A" A" 1 1 

- SFm2m{x2,yi)^{SFi'y5Y)^2rm{y2,Xl)e"'''''' + {h5Y^F)n2mi{y2, Xl)e''''y^ j , 

M'^ISfKz) = -mo J rf^xj^e-^P(^-^V(V){tr[275Y'5^(^'^)] +t^['5^(^''2)^75Y]}- (§5) 
To calculate the decay constant, we define the formula 

f - lim rr.-pi^^SiM2) m) 

- P^i^. i^^"-^^'^'^^^55^_(x,y) ^^^^ 

1 _ ^ Jd^J^JSFKO) , , 

P^Pb dSFnm{x,y) 



corresponding to the exact definition (p6D , or 



U := ^hm^— (P|J5"^[^,^](0)|0). (88) 

Therefore systematically approximated can be obtained for any truncation of Fioopi^'i?]. 
The previous formula Eq.(|79|) coincides with Eq.(^) if the rainbow- ladder approximation is 
employed. 



4.1 Pagels— Stokar Formula 

Pagels and Stokar proposed a useful approximation for the decay constant in the chiral 
hmit in terms of the constituent quark mass function B{q^) of the SD equation. ||2^, ^ The 
Pagels-Stokar formula is 

,2_Nc r.2 , B{-ql)[2B{-ql) + qlB\-ql)\ 

^-"Svr^io ^^^^ [ql + B\-ql)f ^^^^ 

where B{—q%) is defined by the ansatz 

instead of the general form of the SD solution 

B'{x) denotes the derivative of B{x) and qE denotes the Euclidean momentum i.e. q]^ = —q^. 

This formula Eq . (|89|) can be derived from an approximated BS amplitude with the Ward- 
Takahashi identity for the axial- vector current. Therefore when the axial- vector current is 
modified, it should also be modified. In this section, we propose a new formula similar to the 
Pagels-Stokar formula in the effective model. 

The BS amplitude is given by 
1 A° 

Xnm(g; 0) = ^{^75^, SF{q)}nm (92) 

from the Ward-Takahashi identity for the axial- vector current in the chiral limit. But this solu- 
tion is not sufficient to calculate /jr, because is related to the derivative of the BS amplitude 
with respect to the total momentum Pb '■= (V P^, P). The amputated BS amplitude(or BS 
vertex) Xnm{q] Pb) is defined by 

P P 

Xnmiq; Pb) ■■= Spl^,{q + -Y)Xn'm'{q; PB)Spl^,,^{q - -^), (93) 

and in the chiral limit, we obtain from Eq.(p2D 
1 A° 

Xnm{q\ 0) = ^{^75^, Sp^{q)}nm- (94) 



19 



Consider an approximation for the BS amplitude 

P P 

1 Pb A*^ Pq 

J TT 

In this paper we call this the Pagels-Stokar ansatz. This is a very useful formula because it 
gives the approximated BS amplitude for the "pion" in terms only of the solution of the SD 
equation. As this specifies the dependence on the total momentum Pg, one can estimate the 
decay constant without solving the BS equation, for instance using our formula Eq.([7UD or 
Eq.(^7|). The result is given by 

X {^75y( ^^~'-^^^^~'^V + ifhql) - fi-q'-m) + E^iq; P)]] (97) 

where no summation with respect to a is taken. This formula Eq. (|97|) is useful to estimate 
while the validity of the Pagels-Stokar ansatz Eq.(^) is not established. Eq.(|97|) reduces to the 
original Pagels-Stokar formula ( p9D after the Wick rotation when (i) one neglects the effects of 
the cut-off function f{() and the local 4-quark interaction to the axial-vector current, and (ii) 
if one assumes the ansatz (Wi) of the SD solution. 



5 Numerical Results for a Concrete Example 

In this section we show some numerical results as an example of our approach. We here employ 
the effective model introduced by Aoki et al.[p. The interaction is given by the one gluon 
exchange with the running coupling constant which depends on the special set of momenta 
according to the Higashijima-Miransky approximation. In the infrared region the running 
coupling constant is assumed to be constant. We do not elucidate the detail of this model here. 
(See Ref. [p9| .) Our choice of the parameters is Aqcd = SOOMeV, tw = —0.5, to = —3.0 and 
the ultraviolet cut-off paramater Ayy = 2.0GeV. We have solved the SD equation and the BS 
equation for the pion in the Euclidean momentum in the chiral limit. As is pointed out in 
Ref.0, the consistency of the SD and BS equations in the CJT action formulation guarantees 
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that the pion becomes massless in the chiral hmit. We have confirmed this in our numerical 
calculation. 

In the rainbow-ladder approximation, the exact value of the pion decay constant in the 
chiral limit is given by Eq.(^) with P§ = 0, which gives 

U = 72MeV. (98) 

This can be first compared with the naive value f^^ in which E"'{q,p) term coming from the 
non-local interaction is omitted from Eq.(|T9D (See Eq.(^)). The result ^is 

= 122MeV. (99) 

This is more than 70% deviated from the true value. Thus we see that the correction term 
E°'{q,p) is essential. In Ref. [Q, the authors normalize the BS amplitude by the true decay 
constant as in Eg . (p2|) . To calculate the decay constant, however, the naive definition 
Eq.(p7D is employed. As a result, their definition of the decay constant gives 



'UU = 94MeV. (100) 
Next we examine the Pagel-Stokar formula. The original Pagels-Stokar formula Eq.(^9|) gives 
/^^ = 96MeV (101) 
which is close to Eq.( |100| ) as is pointed outf^n Ref. 0. Our new formula Eq.(^^ gives 

f^^ = 73MeV (102) 



which almost coincides with Eq. (pq) . Thus we find that the new formula gives rather good 
prediction, although only the SD equation is to be solved in order to calculate Eq.(0). Thus 
this is the most economical way to estimate the pion decay constant. 



The cut-off regularization is applied for the naive value Eq. (M) . 



The value of f^^ is different from Ref.|^, because of the different choice of the parameters ^if, Aqcd and 

Auv- 
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6 Conclusion 



In this paper, we have discussed dynamical chiral symmetry breaking in effective chiral quark 
models of QCD. Because the effective models may contain loop momentum cut-off as well as 
non-local interactions, the conserved axial-vector current is modified accordingly. Then the 
low-energy constants and relations, such as the pion decay constant, the Gell-Mann-Oakes- 
Renner relation, become very complicated. Nevertheless, using the CJT action formulation, we 
have proved that the combination of the SD and BS equations preserve the chiral symmetry 
and that the BS equation bears the NG pion solution. It is also shown that the pion decay 
constant must be defined according to the modified axial-vector current and that such 
satisfies the PCAC relation. 

We have derived a general formula of the pion decay constant in terms of the quark full 
propagator and the pion BS amplitude. The formula is consistent with the loop expansion of 
the CJT effective action. A numerical analysis given in Sec.^ shows that the consistency of the 
SD and BS equations with the chiral symmetry is essential for the low-energy relations. We 
have proposed a Pagels-Stokar type formula which gives the pion decay constant in terms only 
of the mass function B{q^) of the SD equation. We have found that this formula gives a very 
good approximation and therefore saves computation time. 

Our intention is to apply the general formulation given in the present paper to the study 
of the light qq mesons from the chiral symmetry viewpoint. As we employ realistic effective 
models, beyond the NJL model, the chiral symmetry is not trivially conserved. Thus we need 
a consistent approach of the system of the SD and BS equations. We believe that the present 
formulation gives a consistent view of the dynamical symmetry breaking in the effective model 
analyses of low-energy hadrons. 



In a separate paper [^, we study the pion in the improved ladder model with finite quark 
mass. The introduction of finite quark mass breaks the chiral symmetry explicitly. It is impor- 
tant and interesting to investigate how far the chiral symmetry can be applied. For instance, 
the Gell-Mann-Oakes-Renner relation is proved from the chiral symmetry in the mo — > 
limit. Its deviation at finite mo should be studied. Such a study will give an indication to the 
applicability of the chiral perturbation theory. [B^] 
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The axial U{1) symmetry is known to be broken by the anomaly, which may be caused 
by the instanton|T^ configration of QCD. In effective models for quarks, the U{1)a breaking 
may be represented by the instanton mediated interaction. In our study of a realistic model 



in the flavor ^[/(S)!^^, we employ the Kobayashi-Maskawa-'t Hooft interact ion [^T], |3^, which 
consists of a six-quark vertex. The present general formulation can be easily extended to 
such six-quark interaction with an appropriate momentum dependence. Our formula for the 
pion decay constant { ^7\) is applicable while the proxy of the axial-vector current requires 
additional terms from the new interaction. 
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